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Abstract 

Using small-angle neutron scattering and liquid integral equation theory, we relate the structure 
factor of flexible dendrimers of 4th generation to their average shape. The shape is measured 
as a radial density profile of monomers belonging to a single dendrimer. From that, we derive 
an effective interaction of Gaussian form between pairs of dendrimers and compute the structure 
factor using the hypernetted chain approximation. Excellent agreement with the corresponding 
experimental results is obtained, without the use of adjustable parameters. The present analysis 
thus strongly supports the previous finding that flexible dendrimers of low generation present 
fluctuating structures akin to star polymers. 

PACS numbers: 61.25.Hq, 61.12.-q, 61.20.Gy, 82.70.Dd 
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I. INTRODUCTION 



Dendrimers are synthetic branched macromolecules with defined structure Starting 
from a trifunctional flexible monomer (generation 0), subsequent shells of trifunctional units 
are attached in a well-defined manner. Fig. [I] displays a dendrimer of 4th generation with 
attached end groups at the ends of the units constituting the last generation. Thus, a tree-like 
structure is generated, which can be viewed as an interesting intermediate between colloids 
and polymers: dendrimers of low generations exhibit enough degrees of conformational 
freedom and will thus present fluctuating structures in solution. Hence, these structures 
will greatly resemble star polymers or lightly branched polymers. If, on the other hand, the 
number of generations is increased, steric interactions between the groups located at the 
periphery of the molecule must result. In this case, significant back-folding of these groups 
must occur and these structures will exhibit a more homogeneous segmental density. Hence, 
dendrimers of higher generation will resemble dense colloidal structures. 

The average radial density distribution T(r) of dendrimers is still a matter of debate. 
Here, the question arises, whether dissolved dendrimers exhibit their maximum segment 
density in the center or at the periphery of the molecule. The first theoretical analysis 
of the radial density distribution by Hervet and de Gennes came to the conclusion that 
dendrimers have a dense shell and hence a minimum of T(r) at the center of the molecule 
0. Subsequent theoretical studies, however, showed that dendrimers exhibit their maximum 
segment density at the center of the molecule [[| £|, |], |7|, |J. Only if electrostatic repulsion 
operates between the segments of the dendrimer, a dense shell structure is to be expected 

Small-angle scattering methods, such as small-angle X-ray scattering (SAXS) |§ and 
small-angle neutron scattering (SANS) JTOf are suitable to investigate the radial structure 
of dissolved dendrimers. Up to now, however, SANS- and SAXS-studies of dendrimers in 
solution did not yet come to a final conclusion regarding the average radial density distribu- 



tion T(r). A survey of literature may be found in a recent review (TT|. A recent study of a 
dendrimer of 7th generation concluded that the internal structure of the molecule is rather 
uniform, with the end groups being preferably located at the periphery of the molecule |13| . 
This result is in contradiction to the theoretical studies discussed above 0, |], ^ || 0, 
SANS-studies of dendrimers of lower generation, however, showed that the molecules under 
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consideration have approximately a Gaussian density distribution jL4, plj. The chemical 



structure of the dendrimers investigated in Refs. |H| and |H , differ with regard to the 
number of generations and the nature of the end groups. The results can therefore only 
be compared with caution. In this respect, we also mention the simulation results of Ref. 
H where it was found that the conformation of dendritic molecules drastically depends 
on the generation number, with a soft, Gaussian-like profile resulting for small generations 
and a "collapsed core" with soft tails for larger generations. Nevertheless, it must be con- 
cluded that the question of the overall structure of dissolved dendrimers has not yet found 
a generally accepted answer (see the discussion of this point in Ref. 



Small-angle scattering conducted at different particle concentrations may be useful for 
further elucidation of this problem. Neglecting incoherent contributions, the measured scat- 
tering intensity I(q) (q: magnitude of scattering vector q; q = (Air/X) sin(0/2); wavelength 
of radiation; 9: scattering angle) may be rendered as ||, [K| 



I(q) = p d (p-p m yis(q)S(q). (1) 

In eq. ([!]) above, p& = N^/V is the number density of iV d dissolved dendrimers in the sample 
volume V. The difference p — p m is the contrast of the scatterer towards the solvent, whereby 
p is the average scattering length density of the scatterer and p m that of the solvent. Is(q) is 
the form factor of the object, a quantity directly related to the mass distribution of scattering 
centers within the macromolecular aggregate, as will explained below. Finally, the structure 
factor S(q) accounts for all interparticle correlations arising from the interactions between 
those. 



The effective interaction between dendrimer centers is formally defined as follows |12j : the 
centers of mass of the dendrimers are kept fixed at prescribed positions and a canonical trace 
is carried out over all monomer degrees of freedom. This procedure yields a constrained free 
energy, depending on the instantaneous configuration of the dendrimers' centers of mass. 
Then, the effective interaction energy between the dendrimers is —k^T times the logarithm 
of this partition function, where ks is Boltzmann's constant and T is the absolute temper- 
ature. When averaging over the dendrimers' positions by employing the so-defined effective 
Hamiltonian, the thermodynamics of the system remains invariant. In general, this proce- 
dure generates two-, three- and higher-order interactions between the centers of mass of the 
dendrimers. However, we will follow usual practice and limit our considerations to two-body 
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potentials only, that is, we will make the usual pair-potential approximation, introducing an 
effective pair interaction V e s(R) between the centers of mass of two dendrimers separated 
by the distance R. We comment on the accuracy of the pair potential approximation in 
section III. Since we are dealing with particles which have an overall spherical shape, this 
effective interaction is spherically symmetric and depends only on the magnitude of the 
vector connecting the two centers. 

The structure factor S(q) is directly related to the pair correlation function g(R) and 



hence to the total interaction potential of the solute molecules [16]. Once V e s(R) is known, 
the calculation of g{R) and S(q) follows from the solution of any of the rich variety of so- 
called liquid integral equation theories [p!6| . For infinite dilution S(q) = 1. Hence, from 
an experimental point of view, eq. (|l|) shows that Is{q) ma Y be obtained from SANS- or 
SAXS-data which have been suitably extrapolated to vanishing concentration [|l(| . 

The structure factor S(q) gives direct information on the steepness of the repulsion of 
the dissolved particles: in the case of hard spheres, a pronounced structuring of the solution 
will occur, which leads to a marked maximum of S(q) If the dissolved particles interact 
via soft repulsion, the maximum of S(q) is decreased. Recently, the problem of S(q) of 
star polymers has been addressed [IT], Here, it could be shown that star polymers 

may be looked upon as ultrasoft colloidal particles. The predictions of theory have met with 
gratifying success when compared to SANS-data measured for concentrated solutions of star 
polymers 0, |2D|. 

Up to now, only a few experimental studies have addressed the problem of the structure 



factor of dissolved dendrimers. Ramzi et al. |2T] studied concentrated solutions of dendrimers 
by SANS. In the case of uncharged dendrimers, these workers found a marked depression of 
the peak of S(q) despite the strong decrease of the structure factor at low q. They concluded 
that the absence of a peak of S(q) is related to the softness of interaction, which results from 
the high internal flexibility of the dendritic structures. If charges are added to the segments 
of the dissolved dendrimers, a strong peak of S(q) is generated which is clearly traced back 



to the screened Coulombic interactions of the now charged species (cf. also Ref. [0].) A 
SANS-study by Topp et al. pelf , showed that the dendrimers under consideration exhibited 
only a weak maximum of S(q), in particular at lower generation. At high concentrations the 
evaluation of data according to eq. (1) did not lead to meaningful results, which indicates 
a considerable change of the internal structure of the dendrimers in this regime. The dilute 
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regime has been studied in Refs. ||T4| , [15|| . The data of S(q) have been interpreted in terms 
of a simplified model, treating the dendrimers as effective hard spheres. It is hence evident 
that the problem of S(q) of dissolved dendrimers is not yet understood quantitatively. 

In this work, we wish to apply a simple theory which allows us to calculate the effec- 
tive interaction between dendrimer centers using the measured density profile of an isolated 
dendrimer as a starting point. It turns out that the interaction has a Gaussian form. 
Subsequently, we employ it to calculate theoretically the structure factor S(q) of dissolved 
dendrimers at various concentrations using the hypernetted chain (HNC) approximation. 
The consequences of the internal flexibility on the measured structure factor S(q) will be 
discussed explicitly. Moreover, a quantitative comparison with the experimental data pre- 
sented in Ref. |15j will be given, showing excellent agreement. It will be demonstrated that 
the radial density distribution T(r) derived from scattering experiments may directly be 
used to calculate S(q). 

II. THEORY 

The starting point for the theory is the monomer density profile p(r) around the center 
of mass of an isolated dendrimer. The experiments of Ref. [|14| , |T5|| offer a direct access to 
this quantity, through the so-called shape function T(r). The scattering intensity Is(q) from 
a solution of dendrimers at infinite dilution has been written as (see eq. (18) in Ref. |fi~4|| ): 

h(q)=T 2 (q) } (2) 

meaning that T 2 (q) is the form factor of the dendrimer and hence T(q) the Fourier transform 
of T{r). Hence, as already pointed out in Ref. |L4|], T(r) is a "shape function" that describes 



the way in which p(r) varies in space. This shape function is dimensionless, hence we may 
write: 

p{r)=p T(r), (3) 

where po is a constant having dimensions of density (length -3 ) and which will be specified 
now. Eq. @ above, together with eq. (2) of Ref. [14| imply: 



J drp(r) = p V p , (4) 
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where V p is the measured partial volume of the solute molecule. On the other hand, the 
integral of p(r) has to yield the total number of monomers N within the dendrimer and this 
fact together with eq. ([|) imply: 

p {r) = ^T(r). (5) 

The shape function T(r) is determined by an inverse Fourier transform of T(q), the latter 
being given by eq. (20) of Ref. [14]. Ignoring the 'tail' (aq 2 + bq) exp(—dq 2 ) there, we have: 



T i (q) = V p i exp(-q 2 R i g /3). (6) 



From eq. (||) and after an inverse Fourier transformation, we readily obtain: 

It can be easily seen that J drT(r) = V p , in agreement with eq. (2) of Ref. [14" . 
From eqs. (H) and (|7]) we obtain the monomer density as: 

M = N {*^T <*>{-%%)■ < 8) 

We now assume that two such dendrimers are kept with their centers of mass at a separa- 
tion R apart and wish to calculate the ensuing interaction potential. Let us assume that the 
monomer-monomer interaction potential is given by some function v(t\ — r 2 ) where ri and 
r 2 are the positions of the two monomers. In the mean-field approximation, i.e., ignoring 
the correlations and possible deformations of the dendrimers, the total interaction potential 
Kff(R-) can be approximated by a double integral over the unperturbed density profiles times 



the monomer- monomer interaction kernel [24]. This approximation should hold when the 



dendrimers are not too deeply interpenetrating and reads as: 

Kff(R) =11 dY 1 dY 2 p{r l )p{Y 2 - R)u(n - r 2 ). (9) 



We now make the simplest possible assumption and model the monomer-monomer interac- 
tion by a delta function (contact repulsion): 

v(ti - r 2 ) = VoksTSfa - r 2 ), (10) 

introducing the excluded volume parameter t> . Eqs. (|) and (|10|) then yield: 

Kff(R) = v k B T J drp( r )p( r - R). (11) 
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In this approximation, the interaction is proportional to the convolution of the monomer 
density of a single dendrimer with itself. But each density profile has the form 
iV(a/7r) 3//2 exp(— or 2 ), i.e., N times a normalized Gaussian with width parameter a. It 
is known that the convolution of a normalized Gaussian having width a with itself is again 
a normalized Gaussian, but with width a/2. Hence, eqs. (|8p and ( JTT| ) yield the final result: 

^^(^r^M-iP- <i2) 

The effective potential depends only on the magnitude R — |R| of the center-to-center 
separation, because of rotational symmetry. Note that the interaction potential of eq. ([12]) 
above, is identical to the Flory-Krigbaum effective interaction potential between the centers 
of mass of two polymer chains [p5 |. The fact that we are dealing with dendrimers enters 
into the relation between N and R g . For long, self- avoiding chains, R g oc iV 3 / 5 holds; in the 
present study this is not the case and we determine both quantities from experiment. 
The important feature is that the resulting interaction has a Gaussian form: 

V eS {R)=ee W (-^), (13) 



a 2 



with 



and 



o \ 3/2 

A-(^) v k B T (14) 



\r,- (15) 

The task is now to calculate e from the experimental data. The number of monomers per 
dendrimer, is N = 94, (see Fig. |I|) where the aromatic end-groups are also included in the 
counting as "single monomers" . 

Next we need an estimate for the excluded volume parameter vq. For this, we set vq = /jj, 
where Iq is the monomer length. Taking the realistic value Iq = 0.4 nm for the latter quantity, 
we obtain: 

v = 0.064 nm 3 . (16) 

The above values, together with the experimentally determined gyration radius R g = 1.85 
nm yield: 

e = 10A2k B T. (17) 
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The corresponding value for linear chains is about 2k bT ||27|| . Hence, the energy barrier for 
the dendrimers at hand is about five times higher as that for linear chains, a result that 
is physically reasonable as dendrimers are more compact objects than chains. A general 
study of the structural and thermodynamic properties of a system of particles interacting 
by a Gaussian potential (the "Gaussian core model", GCM), has been presented recently 



. The GCM shows no freezing for e ^5 100 [|26||, hence we conclude that the system of 



dendrimers at hand will remain fluid at all concentrations. 

Next, we express the density of dendrimers in units that are more convenient. Since there 



are dendrimers in the volume V, the density pd of dendrimers is N^/V. In Ref. [141, the 
volume fraction was used which is related to pa by 

<P = p d V p . (18) 

In order to suppress the dependence on V p , it is natural to use as the unit of length the 
parameter a of the pair potential and to introduce a dimensionless measure of the density, 
namely: 

V = JW 3 , (19) 
o 

which may be regarded as an effective volume fraction of the dissolved dendrimers. From 
eqs. (|TJ), (|TJ) and (0) we obtain: 

" = 1(3) if = - 33 ^ < 20 > 



where the known values R g = 1.85 nm and V p = 15.1 nm 3 have been used ||14j| . 

For every value of rj, the Hypernetted Chain (HNC) equation Jl6| for a system interacting 
via (3v(r) = 10.42 exp(— r 2 /cr 2 ) was solved and the structure factor S(q) was obtained as a 
function of qa . Using eq. (|15|) we obtain o = 2.136 nm. This is in close agreement with 
the effective radius of 2.4 nm estimated in Ref. by a simple model valid only for dilute 



solutions. In what is to follow the g-scale was reexpressed in nm 1 units. 



III. RESULTS AND DISCUSSION 



As can be seen from Fig. [|, the Gaussian interaction potential yields excellent agree- 
ment between theory and experiment at all three concentrations measured. Moreover, this 
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is achieved without the use of free fit parameters and it is physically reasonable that soft, 
interpenetrable objects such as dendrimers interact by means of a correspondingly soft in- 
teraction. The boundedness of the effective interaction at zero separations between the 
dendrimer centers is also physically correct: configurations where the centers of mass of two 
dendrimers coincide are possible, without violation of excluded-volume conditions, hence 
the effective potential does not diverge at the origin. This is quite analogous to the case of 
polymer chains, where the effective interaction is known to be Gaussian |27| . 

The precise numerical value e = 10.42/^7" of the energy barrier, eq. ([H]), yields an 
optimal agreement with the experimental data. However, small deviations from this value 
can also be tolerated and by comparing with the SANS data we have established the limits 
in which this parameter can vary as 8.0 < e ^$ 12.0, yielding a corresponding tolerance 
interval 0.370 nm < 1$ < 0.425 nm for the "effective monomer length" Iq. 

A prominent characteristic of the Gaussian potential at the reduced temperatures con- 
sidered here (and also at higher ones), is its property to yield structure factors which do not 
show any pronounced peak with increasing density, in direct contradistinction with hard, 



diverging interactions such as the hard sphere (HS) potential [26]. In order to demonstrate 
this point, we have also attempted to fit the experimental data using a hypothetical HS 
interaction between the dendrimers, with the HS diameter <ths as a fit parameter. The re- 
sults are shown in Fig. |3|. The best fit at the lowest concentration is obtained by the choice 
o"hs = 1.95.R S , which is thereafter kept constant. It can be seen that the fit quality worsens 
with increasing density, as the HS interaction yields a too high peak as well as a structure 
factor at low g-values which lies below the experimental data. However, as we are in the 
very dilute regime, the strong differences between the structure factors produced by these 
two interactions (Gaussian and HS) are not very pronounced. These differences become 
evident only if one looks at higher densities. The result found here is therefore in agreement 
with the findings of Ref . Jl4j . There it has been found that modeling the structure factor in 
terms of a simple hard-sphere ansatz leads to a satisfactory description of the data in the 
region of lowest volume fractions. 

It is therefore evident that a more stringent test of theory can only be achieved by 
SANS-studies conducted at much higher concentrations. Experiments using the G4 and 
G5 dendrimer used in Ref. [jrj|, 15 are under way. Here we present further evidence for 
the validity of the proposed interaction by comparing with and by offering a theoretical 
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explanation for data already existing in the literature |2T], p23| . 

In Fig. 5 of Ref . [^TJ , Ramzi et al. display the experimentally determined structure factor 
for a very wide concentration range of dendrimers, ranging from dilute to above the overlap 
concentration. It is seen that S(q) is deprived of any significant structure; it has a weak peak 
of height ~ 1.1 at intermediate concentrations and thereafter the peak disappears and S(q) 
becomes a monotonic function of q. At the same time, this phenomenon is accompanied 
by a monotonic reduction of the osmotic compressibility of the solution, the latter being 
proportional to the q — > limit of S(q). This is precisely the behavior of S(q) of the 



Gaussian core model [pq| . The latter has a freezing and reentrant melting transition with 
increasing density at energies e > lOO/c^T, meaning that the liquid S(q) has an anomalous 
dependence on the density: the height of its principal maximum grows up to a certain 
density and then it diminishes again. As the temperature is increased, the anomaly in 
S(q) remains, but the height of the principal peak becomes smaller. In order to provide 
a semi-quantitative comparison with the data of Ref. [21], we show in Fig. [| the evolution 
of S(q) of the GCM with density for an energy barrier e = lOfc^T and packing fractions 
r] = 0.01 — 2.00, corresponding to ^ 0.67. A striking similarity with the results of Ref. 



21| can be easily seen. 

Similar conclusions hold when we compare our predictions with the scattering data of 



Topp et al, Figs. 7 and 8 in Ref. ||23||. There, the S(qYs from samples of two different 



dendrimers are displayed; the same anomalous dependence of S(q) on the concentration 
is seen, with the additional feature that the 5"(g)'s from the dendrimers with the larger 
generation number have higher peaks than their counterparts from the smaller dendrimers 
at the same concentration. In our language, a larger generation number, which implies a 
larger monomer number N, corresponds to a higher energy barrier e, see eq. (|l"4l) . This 
automatically causes S(q) to develop stronger peaks for larger generation numbers. 

Of particular interest is the shape of S(q) for very high concentrations which, as can 
be seen in Refs. |21| and |23| is a monotonic function of q. This behavior can be fully 
understood in the framework of the Gaussian interaction potential. It has been recently 
shown that bounded potentials in general, show at high densities a particular 'mean-field 
behavior' |26], EBL ^9j. The direct correlation function c(r) in this limit becomes equal to 
—/3v(r), where v (r) is the interaction potential at hand and (3 = (fc^T) -1 . For the case of 
the Gaussian potential, this implies that the structure factor at the high density limit takes 
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the form 01 El: 



S(q) = ; — F 5T, (21) 

1 + vr 3 /2^ pa 3 exp [_ (Qa/2) 2 ] 



or, more generally for any bounded potential v (r) 



5(g) = 1— - (22) 

l + /3pu(g)' 



where ■y(g) is the Fourier transform of v(r). Eq. ( |21| ) shows that, for the GCM, 5(g) has 
the form of a "smoothed out" step function. It is small at low g's and approaches unity at 
high g's, as v(q) — > there. The same conclusion holds, actually, for any bounded potential 
whose Fourier transform is a monotonically decaying, nonnegative function of q [f28| , as can 
be seen from eq. (|22|). 

The anomalous behavior of the peak heights of the structure factor would be impossible 
if the assumed interaction between dendrimers were diverging at the origin as any power 
law. Indeed, power-law systems are known to undergo a freezing transition, this implying 
that the corresponding S(q) develops stronger and stronger peaks with increasing density 
until, ultimately, freezing occurs when the principal peak reaches the quasi-universal Hansen- 



Verlet value 2.85 [g(J. Hence, effective dendrimer-dendrimer interactions with such kinds 
of divergence at the origin can be immediately ruled out. This argument, however, does 
not rule out all diverging interactions. Indeed, the peak-height anomaly has also been 
observed in the framework of a theoretical treatment of star polymers fL8"l , whose effective 
interaction is diverging logarithmically at the origin. This behavior, which has also been 



experimentally seen ||31|| , is again intimately related to the reentrant melting behavior of 
these systems |32|, [33]. In fact, it has already been suggested |^3| that dendrimers may 
resemble star polymers below the critical arm number f c = 34 at which crystallization 
is marginally possible [^, |3B|, |34J. However, at variance with the Gaussian interaction, 
the star-star effective potential [|19|] does not have a mean-field high density limit and the 



structure factor of stars never becomes a monotonic function of g, even at extremely high 
densities [ lSj , a feature which appears for bounded interactions only. In this respect, there is 



a difference between the effective interaction between the central monomers of two dendrites 
or two stars (which diverges) and that between their centers of mass (which does not). 

On purely theoretical grounds, the underlying assumptions leading to the derivation of 
the logarithmic interaction between star polymers do not hold for dendrimers. Unlike stars, 
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dendrimers do not obey a power-law dependence of the density profile around the center; 
they show no self-similarity, captured in the Daoud-Cotton blob model of the stars |J]J; and 
they are at least one order of magnitude smaller than star polymers, with the implication 
that many of the notions of polymer physics, based on long chains and universality cannot 
be automatically carried over to dendrimers. 

We conclude this section with a remark on the accuracy of the pair potential approxi- 
mation. Though many-body forces between the dendrimers are necessarily present in the 
solution, there are good reasons to believe that their effect can be neglected to a very good 
approximation. This belief is based on the one hand on the corresponding findings on star 
polymers [[36 1. By a direct measurement of the triplet forces, it was established there that 
these have a very small effect on the total force on a star center, and this only shows up at 
concentrations considerably beyond the overlap value. On the other hand, a recent simula- 
tional work by Louis et al. ]27||, reached a similar conclusion for single chains. There, it was 
found that the pair potential yielding the correct thermodynamics of the system remains 
density-independent up to the overlap concentration and displays only a very weak density 
dependence above the latter. This density dependence is, of course, just another way to 
formulate the many-body effects. Once more, they were found to be of minor significance. 
Barring any dramatic alterations in the conformations of single dendrimers upon increasing 
density (such as collapsing), we expect the Gaussian pair potential picture to capture the 
salient characteristics of the behavior of this system at all concentrations. 



IV. SUMMARY AND CONCLUSIONS 



By using the density profile of an isolated dendrimer and the experimentally determined 
characteristics of the macromolecules as input, we derived an effective interaction between 
the dendrimers' centers of mass which is Gaussian in form. By direct comparison with 
experimental data at dilute solutions, we show that this interaction provides an excellent 
description of the measured structure factor. Moreover, it reproduces correctly the as of 
now unexplained trends and features of experimental structure factors at higher densities. A 
direct comparison with SANS data at high concentrations is, evidently, of crucial importance 
for further putting the proposed theory into a strong test. Work along these lines is currently 
in progress. 
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FIG. 1: Structure of dendrimer G4. 
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q[nm ] 



FIG. 2: Experimentally measured (points) and theoretically calculated (curves) structure factors 
S(q) for solutions of dendrimers at three different concentrations. The theoretical curves have been 
obtained by employing a Gaussian effective interaction between the dendrimers. 
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1.1 





FIG. 3: Same as Fig. ^, but now the theoretical curves have been obtained by employing a hypo- 
thetical hard sphere (HS) interaction between the dendrimers, with a HS diameter <ths = l-95i? 9 . 
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0.0 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 
0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 

FIG. 4: The structure factor of the Gaussian core model at energy e = as a function of the 

"packing fraction" rj given in eq. @. From left to right: rj = 0.01, 0.02, 0.05, 0.10, 0.20, 0.40, 
1.00 and 2.00. 
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